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By means of scanning tunneling spectroscopy (STS) we study the electronic structure of the III-V 
semiconductor surface InAs(lll)A in the field emission regime (above the vacuum level). At high 
sample bias voltages (approaching +10 V), a series of well defined resonances are identified as the 
typical Stark shifted image states that are commonly found on metallic surfaces in the form of held 
emission resonances (FER). At lower bias voltages, a more complex situation arises. Up to three 
double peaks are identihed as the hrst three FERs that are split due to their interaction with the 
periodic surface potential. The high corrugation of this potential is also quantihed by means of 
density functional theory (DFT) calculations. Another sharp resonance not belonging to the FER 
series is associated with an unoccupied surface state. 


PACS numbers: 73.20.At, 74.55.+V, 68.37Ef, 68.47Fg 

I. INTRODUCTION 

At any solid-vacuum interface in which a change of 
polarizability takes place, a series of Rydberg-like un¬ 
occupied electronic states exists. They are called im¬ 
age states and result from the coulombic tail of the po¬ 
tential in the vacuum side. Their energies lie below 
the vacuum leveP^ and they are physically located in 
the near-surface region. Their wave functions are con¬ 
fined in the direction perpendicular to the surface, but 
disperse freely in the direction parallel to the surface. 
There are a number of suitable experimental techni^es 
to detect them, like low-energy electron diffractiorP or 
inverse photoemissioiP. By means of scanning tunnel¬ 
ing spectroscopy (STS) these states can also be mapped 
with atomic spatial resolution. The electric field created 
by the STM tip strongly modifies the interface poten¬ 
tial, and the energy at which these states are detected is 
shifted to the field emission regime (above the vacuum 
level) due to the Stark effecfP. 

The existence of these states has been exploited in 
the past for many purposes, like obtaining chemical 
contrasfP, measuring work function fluctuation^, achiev¬ 
ing atomically resolved scanning tunneling microscopy 
(STM) images on insulator^ or even as qubits for 
quantum computatioiP. In recent years, increasing ef¬ 
forts have been made to understand how these elec¬ 
tronic states are modified in the presence of s urface 
nanostructure j^^M^ , steps and surface defect^^^^ or rip¬ 
pled graphenej ^^ * ^^ * to give just a few examples. 

Although much experimental and theoretical research 
has been conducted on image states on metals, compar¬ 
atively few st udies have been performed on surfaces of 
semiconductor^^^I^H^. In the present work, we employ 
STM and STS to study the rich electronic structure in the 
field emission regime of the III-V semiconductor surface 
InAs(lll)A, and identify its image states from their ob¬ 


served field emission resonances (EER), which we analyze 
using a simple one-dimensional model potential defined 
between the surface and the STM tip. 

The InAs(lll)A surface exhibits a 2 x 2 indium- 
vacancy reconstruction which corrugates the ideal 1x1 
surface. Using density functional theory (DFT) we quan¬ 
tify the resulting corrugation of the electrostatic poten¬ 
tial near the surface. A consequence of this corrugation 
is that the lower-order image states, whose electrons are 
confined closer to the surface, appear as split double res¬ 
onances due to their interaction with the periodic poten¬ 
tial of the surface. This feature can be understood in 
terms of scattering of the image state electrons by the 
corrugated surface potential, which prevents them from 
dispersing freely in contrast to their behavior at metal 
surfaces. 


II. EXPERIMENTAL AND THEORETICAL 
METHODS 

The experiments were carried out in an ultra-high vac¬ 
uum chamber equipped with a low-temperature STM 
{Createe GmbH) operated at 5 K. The substrate was 
prepared by growing 20 nm-thick undoped InAs layers 
on top of an InAs(lll)A wafer (from Wafer Teehnology 
Ltd.) by means of molecular beam epitaxy (MBE) moni¬ 
tored by in-situ reflection high energy electron diffraction 
(RHEED). To enable sample transfer to the UHV system 
of the STM apparatus under ambient conditions, a cap¬ 
ping layer of amorphous arsenic was deposited, and later 
desorbed by annealing at 630 K right before inserting the 
crystal into the STM head. We used a standard etched 
polycrystalline tungsten STM tip prepared in vacuum 
by front sputtering and annealing. Once the tip was in 
the STM stage, we performed repeated current pulses at 
sample bias voltages of up to 10 V. This treatment gives 
rise to an agglomeration of indium at the tip apex as ev- 
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ident from the tip behavior in subsequent atom manip¬ 
ulation experiments^. The high voltage treatment was 
usually followed by gentle tip-surface contact to form a 
final atomically sharp tip apex. The differential conduc¬ 
tance signal dljdV (/ denotes the tunneling current and 
V the voltage) was measured with a lock-in amplifier 
using a bias oscillation frequency of 675 Hz and an am¬ 
plitude of 40 mV. All bias voltages are referred to the 
sample with respect to the tip. 

First-principles DFT calculations were used to deter¬ 
mine the equilibrium geometry and near-surface elec¬ 
trostatic potential of the InAs(lll)A-(2 x 2) surface. 
Total energies and forces were calculated within the 
generalized-gradient approximation of Perdew, Burke 
and Ernzerhof functional (PBE]P^ to DET, using 
projector-augmented-wave (PAW) potentials as imple¬ 
mented in VASI^. The plane-wave cutoff for all cal¬ 
culations was 250 eV. The surface calculations were per¬ 
formed in a slab geometry with ten layers of InAs and a 
vacuum region of 20 A. The topmost three atomic layers 
were relaxed until the largest force component on ev¬ 
ery atom was below 0.01 eV/A. The sampling of the 
2x2 surface Brillouin zone was carried out with a 3x3 
Monkhorst-Pack mesh centered at the F point. 


III. RESULTS AND DISCUSSION 


The red curve in Fig.j^a) corresponds to the differen¬ 
tial conductance signal dl/dV measured on InAs(lll)A 
for bias voltages between +2 and +10 V. To prevent cur¬ 
rent overload as the bias voltage is ramped up into the 
field emission regime, the spectrum was acquired at con¬ 
stant current (feedback loop switched on), so that the 
tip height was changing continuously. The tip retraction 
was acquired simultaneously and is represented in green 
in this panel. From +3 to +10 V bias, 15 peaks are clearly 
visible in the spectrum, and our aim is to identify their 
nature. On top of the spectrum, a black curve marks the 
fit of the red curve to 15 lorentzian peaks with a cubic 
polynomial background. The individual fitted peaks are 
also shown in their corresponding locations. 

A first inspection reveals that the identification of all 
peaks is not trivial, especially in the low bias range, 
where the energy separation between consecutive peaks 
does not seem to follow a clear pattern. To analyze this 
spectrum we develop here a simple model to help us in 
discerning whether a particular peak assignment is rea¬ 
sonable or not. The aim is to know where to expect 
the presence of the different FERs on InAs(lll)A along 
an energy resolved spectrum such as the one shown in 
Fig. a). We solve the time-independent Schrodinger 
equation for an electron confined in the one-dimensional 
potential barrier shown as a black curve in Eig.j^b). The 
integration was done numericall 3 @l using the Numerov 
method^. The potential is the sum of three contribu¬ 
tions: 


(z) the surface image potential, 

-1 / e - 1 

47reo 4 \e + 1 

(zz) the tip image potential, 

— 1 / 1 

and 

(zzz) the electrostatic potential between sample and tip. 
Assuming a spherical tip, the value of this poten¬ 
tial along the cylindrical symmetry axis (coordinate 
z, representing the distance from the surface plane 
in the direction to the tip) can be calculated as 
described in Ref. [26] and depends on the effective 
radius Rtip, the bias voltage, the tip-sample sepa¬ 
ration and the difference in work function between 
sample and tip, ^sample - ^tip- 

In expression 0, e is the electron charge, eo is the vac¬ 
uum permitivity and e is the sample permitiviW, which 
has a value of 15.15 for InAP^. In expression Z{V) 
is the experimental tip displacement (green curve in 
Eig.j^a)) with respect to the initial tip-sample distance 
Zq corresponding to a set point of 0.1 V (initial bias in 
the spectrum) and 0.1 nA (constant current during Z{V) 
spectroscopy). Zq has an estimated valu^^of 5.5 A. We 
choose the absolute zero of the potential to be located 
at the vacuum level of the sample. The potential is as¬ 
sumed to be infinite at both sides of the vacuum gap 
(hard-wall boundary conditions), for which we neglect 
the bulk potential in this approximatiorP. This assump¬ 
tion is equivalent to forcing the wave function to be zero 
at both sides of the tip-sample gap. Since image states 
are confined states in the 2 : direction mainly present in 
the tip-surface gap, we expect our model to provide a 
suitable scenario for fitting the experimental data. 

This model has three free parameters, namely the sam¬ 
ple and tip work functions and the radius of the tip. In 
the simulation, the value of these parameters was ad¬ 
justed so that the maximum number of resonances can 
be explained. The vertical lines in Eig.j^a) show the 
fitted EER positions thus obtained using ^sample = 4.3 
eV, (/)tip = 2.8 eV and Rtip = 40 A. The fit reproduces 
well the position of the last eight measured resonances 
and yields reasonable agreement for the first three, which 
we assume to be double resonances (labeled na and n/3, 
where n = 1,2,3). The sharp resonance labeled S does 
not belong to the calculated series of EERs and is dis¬ 
cussed later as well as the origin of the EER splittings. 
Although other peak assignments might be conceivable 
(for example, bulk related states), we discuss only the 
scenario in which these resonances are derived from image 
potential states, provided that a simple potential model 
yields a reasonable agreement with the data. 
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FIG. 1. (Color online) (a) Differential tunnel conductance 
dl/dV (red curve) recorded on bare InAs(lll)A as a func¬ 
tion of the bias voltage and at a constant current of 8 pA. A 
fit of the spectrum to 15 lorentzian peaks plus a polynomial 
background is also shown (black curve). The tip retraction 
with respect to the initial tip-sample distance of about 5.5 
A is shown in green. The blue vertical lines mark the lo¬ 
cation of the field emission resonances (FER) as calculated 
by solving the one-dimensional Schrodinger equation for the 
potential depicted in (b), using the indicated values for the 
parameters Rtip, 0tip and ^sample- (b) One-dimensional po¬ 
tential model for two different tip heights, corresponding to 
the bias at which the tunnel junction is in resonance with the 
FERs n = 1 and n — 2. The associated electron probability 
densities are shown in blue. The corresponding expectation 
value for the electron position along the surface-tip axis {z)^ 
is also indicated. 


The experimental work function of InAs(lll) A has not 
been reported in the literature. Work function values for 
III-V compound semiconductor surfaces typically vary 
between 4.4 and 5.4 eV (Refs.[29l and [30]). For exam¬ 
ple, the reported value for the As-terminated surface 
InAs(lll)B is 4.7 eV (see Ref. [31]). Therefore, the sam¬ 
ple work function of 4.3 eV obtained in the fit appears 
reasonable. On the other hand, the unrealistically small 
tip work funtion of 2.8 eV and the tip radius of ^40 
A should be considered effective values arising from the 
simplifying approximations of the underlying model. 

To improve our understanding of the observed elec¬ 
tronic structure, we modified the local electrostatic po- 
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FIG. 2. (Golor online) (a), (b) Differential conductance maps 
(constant current 10 pA) measured as a function of bias volt¬ 
age and position along the white dashed line marked in the 
corresponding STM image shown above (0.5 V, 0.1 nA). The 
dashed line runs across a single In adatom (a) and a 5- 
membered adatom chain (b). The surface unit cell is marked 
with a white rhombus in the STM maps, c) Line profiles 
extracted from the map in (b) at x = 0 A (in red) and at 
X = 100 A (in blue). The curves are offset along the vertical 
axis for clarity. 


tential above the surface to see its effect on the differ¬ 
ent resonances. On InAs(lll)A, this can be achieved 
by using charged defects such as indium adatoms (In), 
which are natural electron donors with a charge state 
+1 (Ref.pT]). In adatoms (either native or artificially 
created by soft tip in dentat ion) can be readily reposi¬ 
tioned by the STM tipP^^^. Due to the low screening 
provided by the semiconductor surface, the In adatom 
charge is strongly localized and by repositioning and ar¬ 
ranging them we can tailor the local electrostatic poten¬ 
tial. 

The differential conductance maps in Fig. ^a,b) show 
how the different resonances described for the bare sur¬ 
face in Fig.[^ are affected in the presence of two nanos¬ 
tructures: a single adatom and a 5-membered adatom 
chain. The maps detail the bias dependence of dl/dV 
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profiles measured along the dashed lines in the STM im¬ 
ages at the top of Fig.|^a,b) respectively. Due to the 
larger charge associated with the chain, its effect is more 
dramatic. On the basis of this experiment, a number of 
observations can be stated. On the one hand, every dou¬ 
ble resonance measured on the bare surface (labeled with 
numbers 1, 2 and 3) is merged into a single enhanced res¬ 
onance when measured on top of the nanostructure. The 
effect can be seen more clearly in Fig.|^c) where the pro¬ 
files extracted from Fig.j^b) at x = 0 A and at x = 100 
A are displayed for direct comparison. We observe also 
that all resonances are shifting to lower binding ener¬ 
gies, as expected for the lowering of the local potential 
due to the positive charge provided by the nanostruc¬ 
ture. This shift is also visible in the case of resonance 
S. However, contrary to the trend found for the FERs, 
the peak marked with S becomes broader and less in¬ 
tense on top of the nanostructure, confirming its distinct 
nature with respect to the rest of states in the series. 
The electronic structure of InAs(lll)A in the unoccu¬ 
pied region has been studied with inverse photoemission 
spectroscopy (IPES) (see Ref. [33] ). According to that 
study, a surface state exists at ^ 4.6 eV with respect to 
the Eermi level of the surface. This value is fully com¬ 
patible with the binding energy of resonance S obtained 
here with STS. We therefore conclude that S cannot be 
associated with an image state, and should be instead 
assigned to an intrinsic surface state, that broadens in 
the presence of the nanostructure. 

Next, we address the origin of the splitting observed 
for the lower-order EERs. The associated electron prob¬ 
ability densities for the first two EERs is shown in blue in 
Fig-lH These curves represent the calculated electron dis¬ 
tribution along the z direction (perpendicular to the sur¬ 
face) in the sample-tip gap. Note that in the ID model 
used to obtain these distributions, the EER electrons dis¬ 
perse freely in the direction parallel to the surface but 
are confined in the z direction. The expectation value 
for the electron position along the surface-tip axis {z)^^ 
is 2.9 A for the EER with n = 1. This value is quite 
close to the surface, raising the question of whether the 
crystal potential is sufficiently smooth at ^ = 2.9 A to 
ensure free dispersion of the EER wavefunction paral¬ 
lel to the surface. In fact, we find that the potential is 
quite corrugated. Eigurej^a) shows the full DET electro¬ 
static potential 2.9 A above the surface. Eig.|^b) shows a 
cross-section of the potential containing the surface nor¬ 
mal and the in-plane direction at x = 0 marked by the 
vertical dashed line in Eig.j^a). The In and As atoms 
in this plane are labeled. The intrinsic vacancy is lo¬ 
cated at (z^y) = (0,0), where an In atom is missing in 
the topmost atomic layer. The red curve shown below 
the potential map in Eig.|^a) shows the potential along 
the horizontal dashed line. This profile runs across the 
positions of the maximum and minimum in the entire po¬ 
tential map, corresponding to the intrinsic vacancy sites 
and the In surface atoms respectively. A stick-and-ball 
model of the surface bilayer is shown below. 



FIG. 3. (Color online) (a) DFT calculated potential map 
for InAs(lll)A at a plane located 2.9 A above the topmost 
atomic layer of the surface. The black rhombus marks the 
surface unit cell, and each of its vertices is located at the in¬ 
trinsic vacancy sites that give rise to the surface reconstruc¬ 
tion with (2 X 2) periodicity. The red curve shown below is 
the profile across the center of the unit cell, as marked in the 
map with a red dashed line. The black dotted lines mark 
the isocontours for two values of the potential, —0.25 V and 
-h0.15 V for a better visualization of the potential topology 
in the xy plane. Below, the stick-and-ball model of the top¬ 
most bilayer of atoms (indium depicted in green and arsenic 
in red) helps to identify the surface orientation and locations, 
(b) DFT calculated potential map within the x = 0 plane, 
which intersects the map in (a) along the dashed blue line. 
The intrinsic vacancy is located at {z,y) = (0,0), where an 
In atom is missing in the last atomic layer. The isocontours 
for —0.25 V and -h0.15 V are also represented, (c) Potential 
corrugation (Pmax — Pmin) for every plane z, on the vacuum 
side. The potential is still considerably corrugated for rela¬ 
tively large values of z. The blue dashed line in (b) and (c) 
corresponds to the plane z = 2.9 A depicted in (a). 


The corrugation of the DFT potential at any given 2 ; 
plane, which we define as the difference between its max¬ 
imum and minimum (Pmax — Pmin), is very large near 
the surface (0.7 V at 2 ; = 2.9 A) and rapidly decreases 
at larger z, as shown in Fig.j^c). The electronic states 
belonging to the first FERs have a significant presence in 
the near-surface region, where the potential corrugation 
is still high, and thus undergo scattering by the periodic 
potential in the direction parallel to the surface. Due 
to this interaction, the low-order states are energetically 
split. If the interaction is strong enough, the energy dis¬ 
persion of the image state in the plane parallel to the 
surface can deviate significantly from the free-electron 
like dispersion assumed in our ID modeP^. 

This reasoning is supported quantitatively by the data 
in Eig.|^ Panels (a) and (b) show spatial maps of the dif¬ 
ferential conductance measured on the bare surface for 
two different bias voltages, corresponding to the EERs 
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FIG. 4. (Color online) (a), (b) Constant height differential 
conductance maps on InAs(lll)A for the bias voltages corre¬ 
sponding to resonances la and 1/3 shown in Fig.[^a), respec¬ 
tively. The surface unit cell is marked with a black rhombus, 
(c) Line profiles along the dashed lines indicated in panels (a) 
(in green) and (b) (in blue) showing different LDOS distribu¬ 
tions across the surface unit cell. 


la and 1^5, respectively. Each of the four vertices of the 
black rhombus marking the surface unit cell, is located at 
an intrinsic vacancy site, and both Fig.and Fig. [^dis¬ 
play the same surface orientatiorP^. Because the maps of 
Fig.|^a) and (b) were measured at constant tip height, 
they provide a measure of the lateral distribution of elec¬ 
tronic states at the particular energies associated with 
each of the two peaks of the first FER doublet. The 
lower-energy peak la has most of its electronic density 
concentrated at the center of the right half of the unit 
cell (Eig.|^a)). This region coincides with the surface 
location at which the potential deepens, as can be seen 
in Eig.j^a). The isocontour labeled with —0.25 encloses 
the region in which the potential is lower than —0.25 V 
and where the absolute minimum is located. In contrast, 
the state density associated with the higher energy peak 
IP (Eig.^b)) distributes mainly at the center of the left 
half of the unit cell, but extends over the rest of the sur¬ 
face with lower density at the vacancy sites, forming a 
honeycomb structure. In Eig.j^a), the region delimited 
by the isocontours at —0.25 and -f0.15 V extends over 
the same region of the xy plane as the honeycomb struc¬ 


ture. Hence, the probability density of the lower (higher) 
energy state is concentrated in the regions where the po¬ 
tential is lower (higher). This situation resembles the 
simple case of gap formation in the nearly-free-electron 
model, where eigenstates at the Brillouin zone bound¬ 
ary have a high expectation value either between (higher 
energy) or centered around the nuclei positions (lower 
energy) defining the minima of the periodic potential. A 
completely free-electron image state would instead have 
a constant probability density across the xy plane, as 
commonly found on metal surfaces. 

We described in Eig.j^how the double resonances ob¬ 
served on the bare surface become single peaks when 
measured on the adatom nanostructures. In the context 
of our model, this phenomenon can be understood as an 
indication that the potential corrugation is smoother on 
the nanostructures as compared to the bare surface. 


IV. CONCLUSIONS 

We have studied the electronic structure of the 
InAs(lll)A surface in the field emission regime and 
gained a qualitative, as well as semi-quantitative under¬ 
standing of the resonances observed in scanning tunnel¬ 
ing spectroscopy. In contrast to the behavior observed 
on metal surfaces, our results reveal a sizable splitting of 
the lower-order field emission resonances (up to n = 3) 
in the low-bias regime. We interpret this splitting as an 
indication that the lower-order surface image states do 
not disperse freely in the direction parallel to the surface, 
due to the relatively large surface potential corrugation 
predicted by DET calculations. The spatial distribution 
of the electronic states observed at the respective ener¬ 
gies of the split resonances confirms this picture. In ad¬ 
dition to the field-emission-resonance states, the spec¬ 
troscopic data reveal an electronic state that does not 
belong to the series of image states and is assigned to an 
intrinsic surface state of InAs(lll)A. Our results show 
that the surface potential corrugation inherent to a com¬ 
pound semiconductor surface can have significant impact 
on the energy level spectrum of image states in the field- 
emission regime. 
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